Conflict between trajectories and density description: the statistical source of 

disagreement 
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We study an idealized version of intermittent process leading the fluctuations of a stochastic 
dichotomous variable £. It consists of an overdamped and symmetric potential well with a cusp- 
like minimum. The right-hand and left-hand portions of the potential corresponds to £ = W and 
£ = — W, respectively. When the particle reaches this minimum is injected back to a different and 
randomly chosen position, still within the potential well. We build up the corresponding Frobenius- 
Perron equation and we evaluate the correlation function of the stochastic variable £, called &z(t). 
We assign to the potential well a form yielding $|(i) = (T/(t + T)) 13 , with /3 > 0. Thanks to 
the symmetry of the potential, there are no biases, and we limit ourselves to considering correlation 
functions with an even number of times, indicated for concision, by (12), (1234) and, more, in general, 
by (1...2n). The adoption of a formal treatment, based on density, and thus of the operator driving 
the density time evolution, establishes a prescription for the evaluation of the correlation functions, 
yielding (1...2n) = (12). ..{(2n — l)2n). We study the same dynamic problem using trajectories, and 
we establish that the resulting two-time correlation function coincides with that afforded by the 
density picture, as it should. We then study the four-times correlation function and we prove that 
in the non-Poisson case it departs from the density prescription, namely, from (1234) = (12) (34). 
We conclude that this is the main reason why the two pictures yield two different diffusion processes, 
as noticed in an earlier work [M. Bologna, P. Grigolini, B.J. West, Chem. Phys. 284, (1-2) 115-128 
(2002)]. 

PACS numbers: 02.50.Ey,05.20.-y,05.40.Fb,05.70.-a 
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I. INTRODUCTION 



The assumption of the equivalence between trajecto- 
ries and density is one of the major tenets of modern 
physics. In fact, it is of fundamental importance in sta- 
tistical mechanics, where the concept of density was orig- 
inally born, as a consequence of the Gibbs perspective 
resting on the use of infinitely many copies of the same 
dynamic system The concept of density, under the 
form of statistical density matrix is also at the basis of 
the foundation of quantum mechanics and of decoherence 
theory [|| that enables us to interpret all the statistical 
processes occurring in nature as being compatible with 
the unitary time evolution of the whole Universe. As al- 
ready pointed out by Joos years ago ||, the derivation 
of equation of motion for density, non diagonal in the 
position representation of the density matrix, does not 
lead to the unique result of letting trajectories emerge. 
On one side, this seems to make unnecessary to alter the 
quantum-mechanical law of evolution so that localized 
states for macro-objects naturally emerge Q and, on the 
other side, generates the conviction that there is no more 
need for the postulate of wave- function collapse || . 

It comes to be, therefore, a great surprise that the au- 
thors of Ref. p found a conflict between trajectories 
and density picture to exist in a case of non-ordinary 



statistical mechanics. This has to do with the diffusion 
process generated by a fluctuating variable £ (t) , and con- 
sequently with the very simple equation of motion 



dx 
~dt 



£(*)• 



(1) 



The anomalous nature of the resulting diffusion process 
is generated by the fact that the correlation function of 
the stochastic variable is 



T 



t + T 



with 



Q < 0< 1. 



(2) 



(3) 



Notice that the form of Eq. ^ is essential for the realiza- 
tion of anomalous diffusion. For this form of anomalous 
diffusion to be of Levy kind, it is also fundamental to 
assume the variable £ to be dichotomous. In this pa- 
per, we assume, in fact, that £ can only have the values 
W and — W. The Levy nature of the resulting diffusion 
process can be established using either the Continuous 
Time Random Walk § or, as pointed out in Ref. ||, the 
Generalized Central Limit Theorem Q. Both methods, 
however, are based on the direct use of trajectories, and 



2 



consequently are not yet a rigorous derivation of Levy 
statistics from within the Liouville-like approach made 
necessary by the adoption of density picture [jl0|. This 
ambitious task was addressed by the authors of Rcf . jio) , 
but it became clear later fi that the generalized diffusion 
equation built up by these authors becomes compatible 
with Levy statistics through the adoption of a Markov 
approximation that forces it to depart significantly from 
the exact density time evolution. The form of this gen- 
eralized diffusion equation is: 

t) = (£ 2 ) £ - t')-^a(x, t')dt>. (4) 

The main purpose of this paper is to identify the cru- 
cial property that is responsible for the conflict between 
trajectories and density discovered by the authors of Ref. 

However, to reach this important result in a satis- 
factory manner it is necessary to complete the program 
of the authors of Ref. |^|. These authors proposed a 
Frobenius-Perron operator to describe the dynamics of 
the fluctuation process responsible for anomalous diffu- 
sion. We have to prove that this Frobenius-Perron op- 
erator yields the same correlation function as the trajec- 
tories produced by the same dynamics. This is usually 
taken for granted, but in this case, where the equivalence 
between the two pictures is questioned, it is of funda- 
mental importance to prove that density and trajectory 
approach yield the same two-time correlation function. 
This is useful to identify the true source of breakdown. 
As already pointed out by the authors of Ref. |llj (see 
also once the Frobenius-Perron equation is estab- 

lished, we have to work with it without any more direct 
recourse to the concept of trajectories. Thus, we proceed 
as follows. We devote Section 2 to the derivation of the 
Frobenius-Perron operator. In Section 3 we derive the 
correlation function of Eq.(|J) using trajectory. In Sec- 
tion 4 we derive the same correlation function from the 
Frobenius-Perron operator. As earlier stated, this is a 
crucial result, confirming that the source of discrepancy 
between trajectory and density has to be found some- 
where else. In Section 5 we derive Eq.(^J) and in section 
6 we identify the true source of the trajectories-density 
conflict. It is caused by the four-times correlation func- 
tion, and plausibly, of higher order. We devote Section 7 
to concluding remarks. 

II. AN IDEALIZED MODEL OF 
INTERMITTENT DYNAMICS AND ITS 
FROBENIUS-PERRON EQUATION 

As done in Ref. let us consider the following dy- 
namical system. Let us consider a variable y moving 
within the interval / = [0, 2], within a sort of overdamped 
potential with a cusp- like minimum located at y = 1. If 
the initial condition of the particle is y(0) > I, the parti- 
cle moves from the right to the left towards the potential 
minimum. If the initial condition is y(0) < 1, then the 



motion of the particle towards the potential minimum 
takes place from the left to the right. When the parti- 
cle reaches the potential bottom, it is injected back to 
an initial condition in the interval /. This initial condi- 
tion, different from y — 1, is chosen in a random manner 
with uniform probability. We thus realize a mixture of 
randomness and slow deterministic dynamics. The left 
and right portions of the potential V(y) correspond to 
the laminar regions of turbulent dynamics, while ran- 
domness is concentrated at y = 1 . In other words, this is 
an idealization of the map used by Zumofen and Klafter 
, which does not affect the long-time dynamics of the 
process, yielding only the benefit of a neat distinction 
between random and deterministic dynamics. 

The form of the Frobenius-Perron operator is the fol- 
lowing 

§- t p(y,t) = -^(yp(y,t)) + C(t), (5) 

where, as we shall see, 

C(t) oc p(l,t), (6) 

which implies the random injection into any point y of 
the interval /, with equal probability, due to the fact that 
C(t) is independent of y. Eq.(||) must fulfill the following 
physical conditions 

^- [ p(y,t)dy= [ ^-p(y,t)dy = 0. (7) 

at Jl=[0,2] Jl=[0,2] 01 

The condition (f7j) means that the evolution operator pre- 
serves the mass of the distribution that without loss of 
generality is assumed to be given by 

/ p(y,t)dy=l. (8) 

J [0,2] 

To ensure that plugging Eq. (|5|) into Eq. (^) yields the nat- 
ural condition of Eq. (|6|) and, at the same time, the prob- 
ability conservation of Eq. (||) Q , we are forced to adopt 
the following condition for y, 

y = \[Q(l- y )y*-e(y-l)(2-yy}+^Ls(y-l). (9) 

Trandom 

The function 0(x) is the ordinary Heaviside step func- 
tion, A y (t) is a random function of time that can get any 
value of the interval [— 1,-1-1], and T ra ndom is the injec- 
tion time that must fulfill the condition of being infinitely 
smaller than the time of sojourn in one of the laminar 
phases. Eq.(^|), as well as Eq.(H), is the combination of 
two processes, one deterministic and regular, the motion 
within the laminar phase, and the other, expressed by 
the second term on the right hand side of Eq.(^), totally 
random. 

For the calculations that we shall make in Section 4, it 
is convenient to split the density p{y, t) into symmetric 
and antisymmetric part with respect to y = I, 

P(y,t) =p s (y,t)+p A (y,t). (10) 
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Of course we always have that: 

p A (y,t)dy = 



[0,2] 



[0,2] 



Ps{y,t)dy = 1. 



(11) 



The Frobenius-Perron equation of Eq.(|^), with y given 
by Eq. (j^) , expressed in terms of the symmetric and anti- 
symmetric part, yields the following two equations 



0_ 

di 



ps{y,t) 



-A6(l 

+xe( y 



y)g-[y z ps(y,t)} 



y) z Ps(y,t)] + C(t) 
(12) 



and 



0_ 

di 



PA(y,t) 



-A6(l 



-A6(y 



y)g^[y z PA(y,t)} 



y)'p A (v,t)]. (13) 



Note that Eq.(12) yields as an equilibrium state the 
invariant distribution po(y), whose explicit expression is 
easily obtained from Eq.(12) to be 



Po(y) 



2 - z 



6(1 -y) , Q(y-1) 



y~ 



(2-y) 2 



(14) 



The reader can easily check this to be the solution of 
Eq.(|l2|), by setting C{t) = Xp {l), Ps(y,t) = po(y), and 
using Eq.(0) to determine Po(l) and Po{y)- 

We denote this equilibrium state with the quantum-like 
symbol \po). We refer ourselves to a notation where the 
corresponding left state, (fio\ is nothing but the constant 
1. We write the correlation function "^(i) under the 
quantum-like form 

W 2 $s(t) = (i(t)£(0)) = (po||exp(ft)||p ), (15) 

where W 2 = (££)• We note that when equilibrium is 
reached, Eq. (|l2|) leaves room only for the state \po). All 
the other excited states would imply an out of equilib- 
rium condition. Note that the diffusion process under 
study in this paper is stationary. This means that the 
correlation function is evaluated in the correspondence 
of the equilibrium condition \po). The consequence of 
this is significant. In fact the application of the opera- 
tor £ to \po) generates an antisymmetric state. The time 
evolution of this excited state is driven by the operator 
f defined by Eq.(||). This means that the time evolu- 
tion takes place in the antisymmetric space. When we 
apply to this time dependent state the operator £ we 
obtain again a symmetric state. However, the only sym- 
metric state available is (po\. This explains the form of 
Eq. (|l5|) . If we apply the same argument to the four-times 



correlation function we obtain, in the time ordered case 
h < ta < h < k, 

(^1)^2)^3)^4)) = (i(t2-h)i(o)){i(u-t 3 )i(o)}. 

(16) 

This is the four-times case of a more general property 
that, using a concise notation, with an evident meaning, 
we write as 



(12. ...2n) = (12)...(2(n- l)2n). 



(17) 



The authors of Ref. [[| proved that the property of 
17|) yields the generalized diffusion equation of 
. This confirms that Eq.(|^) is a legitimate and 
rigorous consequence of the Frobenius-Perron picture. 



III. THE CORRELATION FUNCTION OF THE 
DICHOTOMOUS FLUCTUATION FROM THE 
TRAJECTORY PICTURE 

Let us focus our attention on Eq.(||) and consider the 
initial condition j/o £ [0, 1]. Then, as it is straightforward 
to prove, the solution for y < 1 is: 



y(t) =yo(l- X(z - l)^" 1 *) 



z _ 1 -l/(z-l) 



(18) 



Using Eq.( 
at which t 
solution is: 



18|) and setting y(T) = 1, we find the time 
le trajectory reaches the point y = 1. The 



T = T(y ) 



1-Vo 



A(* - i)y z - 



(19) 



Note that the function is determined by £(t) — 
i(y(t)). Thus, from Eq.© we obtain: 



w 



[Q(l-y )e(T(y )-t) 



-e(y -l)Q(T(2-y )-t)] 



^2 sign 



\k=0 J 



\fc=0 



\k=0 



(20) 



where the times of sojourn in the laminar phases Tj's read 



ro=T(y o ) = T7 L^o__ 0(1 _ yo) 



1 - (2 - y Q y- 1 
X{z-l){2-y y- 

l-[l + A y ( n )} 



T e( yo -i) 



)pr©(- A ,(n)) 



X(Z-l)[l + Ay(Ti) 

1 - [l - A y (n)Y-i 
+ A(z-l)[l-A y (^- 1 ° (A ^ ri)) - 



(21) 
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Then, we have: 



mm 

w 2 



([6(1- y )e(T(y ) - t) 



+Q(y -l)e(T(2-y )-t)}) + 



+ OC 



E ( si S n (yo - l)sign 



E 

\k=0 



Tk 



'i+l 



e 2>-t -e 5>-t 



U'=0 



\fc=0 



(22) 



As pointed out in Section 2, the calculation of the cor- 
relation function rests on averaging on the invariant dis- 
tribution of Eq. (|lj) . As a consequence of this averaging 
the second term in(p2|) vanishes. In fact, the quantity to 
average is antisymmetric, whereas the statistical weight 
is symmetric. 

It is possible to write the surviving term as: 



mm 

w 2 



l-y* 



2-z) I — fr-— r-t I -— 



1 



with 



/■(l+A^-l)*)- 1 /'*- 1 ) 
- (2 - z) / y~ z+1 dy 
Jo 

= (1 + A(z - l)*) - * 2- *^ -13 
= (l + X(z-l)t)- , 



z - 1 



(23) 



(24) 



Since we focus our attention on < /3 < 1, we have to 
consider 3/2 < z < 2. Note that the region 1 < z < 3/2 
does not produce evident signs of deviation from ordinary 
statistics. However, as we shall see in Section 5, an exact 
agreement between density and trajectory is recovered 
only at z = 1, when the correlation function becomes 
identical to the exponential function exp(— At). Note also 
that Eq.(^3|) becomes identical to Eq.(||) after setting the 
condition 



X(z-l) 



1 

T' 



(25) 



IV. THE CORRELATION FUNCTION OF THE 
DICHOTOMOUS FLUCTUATION FROM THE 
DENSITY PICTURE 

The result of the preceding section is reassuring, since 
it proves that the intermittent model we are using gener- 
ates the wanted inverse power law form for the correlation 
function of the dichotomous variable £(t). However, it is 
based on the adoption of trajectories. In this Section we 
get an even more important result, this being the fact 
that the same correlation function is derived from the 
use of the Frobenius-Perron equation of Eq. (|^) , with no 
direct use of trajectories. 



To fix the ideas, let us consider the following system: 
a particle in the interval [0, 1] moves towards y — 1 fol- 
lowing the prescription y = Xy z and when it reaches 
y = 1 it is injected backwards at a random position in 
the interval. The evolution equation obeyed by the den- 
sities defined on this interval is the same as Eq.(^J), with 
C(t) = Ap(l,t). This dynamic problem was already ad- 
dressed in Refs. 11, [13), and solved using the method 
of characteristics detailed in Ref. |14"| . Let us remind 
the reader that the solution afforded by the method of 
characteristics, in this case yields: 



1 - 



,o 



i 



gziy'-H)- 



where 



5z ( 2 :) = (l + A(z-l) a ;) z /( z - 1 ). 



(26) 



(27) 



To find the quantity (£(t)£(0)) using only density, we 
have to solve Eqs.(|l2|) and (|13|), which are equations of 



the same form as that yielding Eq.(26). For this reason 
we adopt the method of characteristics again. To do the 
calculation in this case, it is convenient to adopt a frame 
symmetric respect to y — 1. Then, let us define Y = y—1. 
Using the new variable and Eq. (p6|) , we find for Eqs.(|H 
and (|l^) the following solution 



Ps(Y,t) 



Ps{0,0 



o «?»((*- 0(1 -l^l)*- 1 ) 
l-\Y 



Ps 



1 



g,((l-\Y\)>-H) 
1 



x 



g z ((l-\Y\r-H) 



(28) 



and 

PA(Y,t) = ii \ [ si»n(y) 
1 



1 - 



l-\Y\ 



g,((l-\Y\>-H)) 



g z ((l-\Y\)*-H) 
Then, the solution consists of two terms: 



,0) x 
(29) 



• the former is an even term and is responsible for 
the long-time limit of the distribution evolution; 

• the latter is an odd term and it disappears in the 
long-time limit. Let us point out why this feature 
is beneficial: Eq.(p"3|) does not contain the injection 
term C(t) and the equilibrium distribution func- 
tion is given by the density distribution of Eq. ([l4]) , 
which is an even function, no matter what the sym- 
metry of the initial distribution function is. 
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As pointed out in Sections 2 (see Eq.(|15|) and Eq.([16|)), 
the correlation function of the operator £ is determined 
by the time independent operator driving the time evo- 
lution of the anti-symmetric distribution. As proved in 
Section 5, this property is of fundamental importance for 
a rigorous derivation of the generalized diffusion equa- 
tion of Eq.(||). In fact, this property means that the 
anti-symmetric space is the irrelevant space, to be traced 
out for the derivation of this important diffusion equa- 
tion. The direct calculation of this correlation function 
yields: 



(|(t)£(0)> = (2-z) 



1 



1 -dY 



£=i- 



gzUi-Y)*- 1 *) 



(30) 



g,((l-Y)>-H) 

The integral (^0|) is exactly solvable and leads to: 

(i(t)i(0)) = (1 + X(z - i) t )-(*-*)/(*-i) (31) 

which is the same result as that found in Section 3, using 
the trajectories. In a similar fashion, it is possible to 
calculate (Y(t)Y(0)): its temporal behavior is a power 
law with the same exponent as that of Eq. (^Tj) . 



V. THE GENERALIZED DIFFUSION 
EQUATION REVISITED 

In this section we shall derive again the generalized 
diffusion equation of Eq. (Q) in a way that is more closely 
connected with the formalism of Sections 2 to 4. We 
define 



\Pi) 



W 



bo). 



(32) 



This is an antisymmetric state. Consequently its time 
evolution is driven by the operator T defined by Eq.(13). 
We now assume |H| that the operator f has a complete 
basis set of eigenvectors with a non-necessarily discrete 
spectrum, even if adopt for simplicity a discrete nota- 
tion, setting aside the easy problem of moving from the 
discrete to the continuous condition. Let us assume 



and 



(A|r+ = (^;. 

so that the biorthonormal condition 



(33) 
(34) 
(35) 



applies. All this allows us to write the correlation func- 
tion "^(i) under the following attractive form 



= (Pi|exp(f*)|pi) =5^(pi|/i)(A|pi>exp(7/.*) 



We are now in a position to derive the generalized 
diffusion equation of Eq.(Q) from a Liouville-like equa- 
tion. It is evident that the Liouville-like equation of 
the whole system, corresponding to the Frobenius-Perron 
treatment of the system of variables £ and y of the pre- 
ceding sections, is given by 



d_ 

df 



H(a:,f,y,t) 



^+r) rn.r.i. ,/./). c57i 



Let us adopt the quantum-like formalism and let us set 

p» = (A|n(t)> (38) 

and 



a(x,t) = p (x,t) = (p |n(i)}. 



(39) 



Thus, by multiplying Eq.(p7|) by (po| and (Al; we get 

d d 
-a(x,t) = -Wj2(pi\rifo PmOM) (40) 

and 



d d 

— Pn(x, t) = -W—a(x, t) + 7 M p M (x, t), (41) 
respectively. The solution of Eq. ([y]) is 

Pll (x,t) = -W J exp(^(t-t'))^-a(x,t')dt'. (42) 

By plugging Eq.(^2|) into Eq.(^0|) we derive the funda- 
mental equation of Eq.([|). It is easy to prove that 
this equation yields a hierarchy of moments (x 2n (t)), with 
n = 1,2, ...,oo, which coincides with the hierarchy that 
would be generated by a fluctuation £ (t) with the correla- 
tion functions obeying the prescription of Eq.([l7|). In the 
next section we shall see that this conflicts with the pre- 
scriptions generated by the adoption of trajectories. The 
arguments of Section 6 are limited to the case of four- 
times correlation function. However, we think that this 
is enough to identify the source of the conflict between 
trajectories and density. 



VI. FOUR-TIMES CORRELATION FUNCTION 

In this section we show that the four-times correlation 
function stemming from probabilistic arguments and tra- 
jectories do not fulfill Eq.(16). We shall see that we can 
write the four-times correlation function in terms of a dif- 
ferent combination of two-times correlation function. In 
this way we can connect higher-order correlation function 
to the waiting time distribution tp(t) previously defined. 
Note that the waiting time distribution and the correla- 
tion function are linked via [fL6| 



(36) 



OO 

ta-h 



(43) 
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This means that the correlation function can be viewed 
as the probability of finding t\ and t 2 in the same lami- 
nar region. In fact, the occurrence of one or more tran- 
sitions, between t\ and t 2 , would make the correlation 
function vanish due to the averaging on uncorrelated fluc- 
tuations, either positive or negative. Using similar argu- 
ments the four-times correlation function can be viewed 
as the probability that both t\ and t 2 belong to one lam- 
inar region and so do £3 and £4. In the following we shall 
use the notation p(ij) as the probability that ti and tj 
belong to the same laminar region. We shall also use 
the notation p(ij) = 1 — p(ij) as the probability that 
at least one transition occurs between ti and tj. With 
the usual Bayesian notation we indicate with p(A, B) 
the joint probability of events A and B and with p(A\B) 
the conditional probability of A given B (thus implying 
p(A, B) = p(B)p(A\B)). It is easy to cast the four-times 
correlation function in our notation as (with W 2 = 1): 



<£(*i)£(*a)£(*s)£(*4)>=p(12,34) 



(44) 



: p(14) + p(23)p(12|23))p(34|23) 



= p(14) + 



(p(12)-p(12,23))(p(34)-p(34,23)) 
l-p(23) 



$ 5 (*4 - *0 



(^(t 2 - tj) - $ c (t 3 - h))(^(u - h) - ^(u - h)) 

where we used the general property: 

P (A)-p(A,B) 



p(A\B) = 



l-p(B) 
and the fact that in our problem 

P(ij,jk) =p(ik). 



(45) 



(46) 



The leading order of Eq. (^4|) is the first term, namely 
<&{(t4 — t\). This yields, for the fourth moment 

t ti ts £2 

(x 4 (t)) = 8 J dt 4 J dt 3 J dt 2 J cM^iKfeKfeK^)}, 



(47) 

and, if the correlation function has a slow decay as in 
(|), (x 4 (t)} oc t 4 -?. Notice that if condition @ held 
true, we would have (x 4 (t)) oc t 4_2/3 . Generalizing the 
4-times treatment to the 2n-times, the trajectory pre- 
scription yields 



(x 2n (t)} oc t 



2n-0 



(48) 



as recently shown in Ref. [|17|. On the other hand, 
the factorization of the even correlation functions, which 



is exact in the density treatment, leads to (x 2n (t)) oc 
t 2n(i-^/2) j f u ifm ing (g), Tms i atter property reflects the 

fact that the solution of the generalized diffusion equa- 
tion (|J) has a solution which rescales with time, while its 
trajectory counterpart does not obey an exact rescaling 

Remarkably, in the exponential case, with some tedious 
algebra it is straightforward to show that 

(£(*l)£(*2)£(*3)£(t4)> = <^l)£fe))(£fe)^4)>, (49) 

which coincides with prescription (|l6|). With boring but 
straightforward arguments it is possible to prove that 
Eq.([49"|) yields the more general condition of Eq.(p"7|). In 
conclusion, in the exponential case the agreement be- 
tween density and trajectories is complete. 



VII. CONCLUDING REMARKS 

The discrepancy between trajectories and density is a 
surprising result. The first report of this disconcerting 
property was given in Ref. ||. However, in that pa- 
per the derivation of the generalized diffusion equation, 
Eq.(^), from a Frobenius- Perron picture was not as ac- 
curate as in the present paper, thereby leaving room to 
the criticism that this equation is not compatible with 
a rigorous Frobenius-Perron picture. This is the reason 
why Sections 2, 3 and 4 have been devoted to a careful 
discussion of the Frobenius-Perron picture. Of remark- 
able interest is the content of Section 4 which derives the 
crucial correlation function of Eq.(^) without using the 
concept of trajectory. There is no doubt, therefore, that 
the density picture adopted is the proper representation 
of the process under study. 

The importance of the Gibbs picture for statistical me- 
chanics is well known, and it is well known that the er- 
godic properties have been introduced to replace the av- 
erage over infinitely many copies of the same system with 
a more realistic average in time iQ. However, there is 
no guarantee that the two pictures are equivalent, in gen- 
eral. Thus, to a first sight the reader might think that 
the breakdown of the equivalence between trajectories 
and density rests on the failure of the ergodic condition, 
a relatively trivial property. The effect discovered by the 
authors of Ref. || does not have anything to do with 
the breakdown of the ergodic condition. If we focus our 
attention on the variable y and on the Frobenius-Perron 
picture of Section 2, we do not find any disagreement 
between trajectories and density, the results of Sections 
3 and 4 lead us to expect that the numerical study of 
a bunch of trajectory yields for the time evolution of a 
non equilibrium probability distribution the same time 
evolution as that prescribed by the Frobenius-Perron op- 
erator . On top of that the variable y is ergodic and the 
generalized diffusion equation of Eq. (||) is based on the 
correct assumption that the Gibbs averages are equal to 
the time average on time of a single trajectory. This is a 
correct assumption. The reason of the failure of Eq. (0) 



7 



to recover the results of numerical calculations based on 
trajectories is much deeper and subtler than the failure of 
the ergodic condition. This reason has to do with the fact 
that the higher-order correlation functions are expressed 
in terms of the two-time correlation functions in two dis- 
tinct ways, the trajectory and the density way. We have 
shown in details that the four-times correlation function 
compatible with the generalized diffusion equation of Eq. 

is the product of two second-order (two-times) corre- 
lation functions. In the non-exponential case, the four- 
times correlation function generated by the trajectories 
breaks this condition. Extending the kind of arguments 
used in Section 6 to correlation functions of higher or- 
der involves much more complicated calculations, but the 
ballistic contribution always shows up, and the formulas 
can be expressed in terms of pairs of two-time correlation 
functions, although by means of expressions of increasing 
length. In conclusion, even if for simplicity it is not done 
here, it is possible to prove that in the case of exponential 
relaxation, the factorization condition is fulfilled by the 
correlation functions of any order. This means that the 
exponential case is the only one where the equivalence 
between trajectory and density is complete. 

What might be the consequence of this result? In the 
introduction we mentioned the possible consequences on 
the foundation of quantum mechanics. One might be 
tempted to conclude that the foundation of Levy pro- 
cesses on the basis of classical trajectories rules out the 
possibility of deriving Levy diffusion from a quantum 
mechanical treatment, in conflict with the conclusions 
of several papers that suggest the opposite to be true 
[||, |o[ §||. We have to point out that all these 
papers studied the kicked rotor, and with this dynami- 
cal system the condition of ballistic motion for extended 
times is realized by a trajectory sticking to the border be- 
tween chaotic sea and acceleration islands. This sticking 
is actually due to the fact that the classical trajectories 
keeps exploring fractal regions of decreasing size. The ar- 
guments used to justify the birth of anomalous diffusion 
rest on classical trajectories, and the quantum breakdown 
of this condition of anomalous diffusion is due to the 
broadening of the quantum wave function that prevents 
it from exploring border regions of even infinitesimally 
small size. It is plausible that the realization of Levy 
statistics cannot be complete, and that the breakdown 
of anomalous diffusion takes place in the very moment 
when tunneling processes are activated. All this seems 
to prove that the papers of Refs. |l9|, |2(], |2l], ^2) cannot 
be used to argue against the conclusion of this paper, and 
cannot be used to support them either. 

We think that the conclusions of this paper might be 
of interest for quantum mechanics at a different level. 
The experiment on resonant fluorescence on single atoms 
have revealed experimentally the occurrence of quantum 
jumps Nevertheless the discovery that the Bloch 

equations consists of two contributions, one compatible 
with the unitary nature of quantum mechanics, and one 
equivalent to the von Neumann measurement postulate, 



makes it plausible to reach the conclusion that, as a rel- 
evant example, the reader can find in the introduction of 
the book of Ref . Q . This is that the hypothesis of spon- 
taneous collapses is unnecessary. The contraction on the 
irrelevant degrees of freedom would be enough to mimic 
properties that are naively interpreted as a manifestation 
of collapses. 

This view is generally accepted, but it is challenged 
by new discoveries. We note that the unraveling process 
compatible with the occurrence of quantum jumps rests 
its foundation on the connection between random walker 
trajectories and master equation, established years ago 
by the authors of Ref. |25| . These authors proved that an 
exponential waiting time distribution implies a Marko- 
vian master equation. Using the more modern jargon 
of de-coherence literature |2j| we can say that the ex- 
ponential waiting time distribution is compatible with 
the quantum mechanical master equations that have the 
Lindblad structure. However, recent experimental stud- 
ies on semiconductor quantum dots p6| , p7| reveal in- 
termittent signals corresponding to a bright state emit- 
ting many photons and to a dark state with no photon 
emission. The waiting time distributions of these two 
states turn out to exhibit inverse power law properties for 
several decades. There are already conjectures that this 
strong deviation from an exponential behavior might be 
due to the existence of slow modulation of barrier widths 
or heights j2§. 

It is worth noticing that the authors of Ref. [^9| used 
this modulation hypothesis to discuss an ideal physical 
problem that it is very similar to the blinking semicon- 
ductors. These authors studied the case where a given 
Pauli matrix a x is forced by the environment to yield a 
symbolic sequence of numbers, corresponding to its two 
eigenvalues, a sequence that can be thought of as cor- 
responding to that of a blinking semiconductor, with, 
for instance, the eigenvalue 1 standing for light, and the 
eigenvalue —1 for darkness. This fluctuating dipole drives 
the motion of another 1/2 spin system, serving a detec- 
tion purpose. Actually, this second spin system is essen- 
tially equivalent to the Kubo stochastic dipole j30j used 
by the authors of Ref. |§. 

The authors of Ref. [^9| found a surprising effect. 
They studied the relaxation properties of the detecting 
dipole driven by the intermittent signal, with a modula- 
tion assumption about the origin of its non-Poisson na- 
ture that fits the conjectures made to explain the prop- 
erties of blinking quantum dots [^8). The waiting time 
distribution was assumed to have a diverging second mo- 
ment, aproperty shared indeed by the blinking quantum 
dots P6| . The calculation of these spectroscopy prop- 
erties was done in two distinct ways. The former was 
compatible with the density prescription, with the only 
constraint that the fluctuation correlation function has 
to coincide with that of the symbolic sequence realized 
by the environmental induced process. The latter cal- 
culation was done, on the contrary, noticing that the 
spectroscopy property can also be expressed through a 
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characteristic function. The explicit analytical form of 
this characteristic function is determined by the central 
limit theorem, either normal, in the case of ordinary sta- 
tistical mechanics, or generalized, in the case of interest, 
that is, the case of waiting time distributions with a di- 
verging second moment. This latter way of doing the 
calculation yielded a result conflicting with the former, 
with the detecting dipole correlation function decaying 
as an exponential rather than as an inverse power law. 
We think that the present article affords a satisfactory 
explanation of why the two methods adopted yield dif- 
ferent results. This is so because only the exponential 
relaxation make the correlation function stemming from 



the density picture identical to those produced by trajec- 
tories. A deviation from the Poissonian condition seems 
to have dramatic consequences. We hope that this paper 
might serve the purpose of triggering further investiga- 
tion on this important issue. 
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